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Abstract 

We show that the Lagrange interpolation polynomials are biorthogonal with respect 
to a set of rational functions whose poles coinicde with interpolation points 

o 

The Newton-Lagrange interpolation is a well-known problem in elementary calculus. Recall 
^ ' basic facts concerning this problem [6], [2]. 

Let A k , k = 0, 1, 2, . . . and a k , k = 0, 1, 2, ... be two arbitrary sequences of complex numbers 
(we assume that all a k are distinct a k ^ aj if k ^ j. By interpolation polynomial we mean a 
Q\ ! n-degree polynomial P n (z) whose values at points a , at, . . . , a n coincide with A , Ai, . . . , A n , 

i.e. 

<; P n (a k ) = A k , fc = 0,l,2,...,n (1) 
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Usually the parameters A k are interpreted as values of some function F(z) at fixed points a k , 



i.e. 



A k = F{a k ) (2) 

In this case polynomials P n (z) interpolate the function F(z) at points a k . 

Explicit expression for interpolation polynomial P n (z) can be presented in two forms. In 
the Newtonian form we have [6], [2] 

00 

p n{z) = ^2[a ,a 1 ,...,a k )uj k (x), (3) 

k=0 

■ 

O where 

^ ■ w = 1, u k (x) = (x- a )(x - ai) . . . (x - a fe _i) 



and [ao, ai, . . . , a k ] denotes the k-th Newtonian divided differnce which is defined as 
[a ] = A , [ao, a x ] = — — — , . . . , [a , a Xl . . . , a k ] = X] — A " 



rV ai - a s=0 u k+1 {a s ) 

5h ' 



where 

k 

w 'k+i( a s) = (a s - ai)(a s - a 2 ) . . . (a s - a s _i)(a s - a s+1 ) . . . (a s - a k ) = J| (a s - aA 

If representation (2) holds where F(z) is a meromorphic function then the Hermite formula 
is useful [2] 

,,,_!/■ F(()d( 
[a ,a 1 , . . . ,a k \ = (2m) / — -, (4) 
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where the contour T in complex plane is chosen such that points ao,ai, . . . , lie inside the 
contour whereas all singularities of the function F(z) lie outside the contour. 
In the Lagrangean form we have [6], [2] 

Introduce the monic interpolation polynomials P n (z) = P n (z)/a n , where 

otn — [ a o? Oi, . . . , a n ] 

In what follows we will assume that a n ^ for all n = 1,2,.... This condition guarantees that 
polynomials P n (z) are indeed of the n-th degree. It is easily seen that P n (z) = z n + 0(z n ~ r ). 
For polynomials P n (z) one has the recurrence relation [3] 

P n +i(z) = {z-a n + -^-)P n (z) - ^-^(z - a n )P n -i(z) (6) 
a n +i a n 

with the initial conditions 

P_! = 0, P (z) = 1 (7) 

It is clear that the set of monic interpolation polynomials P , P\(z), . . . , P n (z), . . . does not 
belong to a set of orthogonal polynomials (OP), becuase OP satisfy 3-term recurrence realtions 
of the form [7] 

P n+1 (z) + b n P n (z) + u n P n ^{z) = zP n (z) (8) 

which doesn't have the form (6). 

Nevertheless, recurrence relation (6) belongs to the so-called -Rj-type recurrence relations 
(in terminology of [4]). It was shown in [4] that polynomials satisfying Rj-type relations possess 
some orthogonality property. In our case this orthogonality property is well known [6]: 

Lemma 1 Polynomials P n {z) satisfy formal orthogonality relation 

I* = i^)' 1 L f = ^ J = 0,l,...,n (9) 
Jr u n+1 (QF(Q a n 

where the contour T encompasses points a , a±, . . . , ajv with N > n and all singularities of the 
function 1/F(z) lie outside the contour. 

For the proof it is sufficient to note, that under conditions upon choice of the contour T, the 
intergral can be presented as a sum of residues 

n a j P (a ) n a j 

I nj = J2 =E rV^' ( 10 ) 

s=0 A s uj n+1 (a s ) s=0 « n w n+ i((i s ) 

where we used interpolaiton property (1). Hence, in the integral (9) one can replace — l/a n 
and we have 

I nj = (2my l / ^ \ = 0, j = 0, 1, . . . ,n - 1 



because the value of the integral from pure rational function doesnot depend on choice of the 
contour T (provided that all poles of the function lie inside the contour) and we can choose T 
as a circle of a great radius. For j = n we have analogously 

I nn = {2-irt) / — = (2m) / — - = 1 a n . 

Jr a n u n+1 (Q Jr a n C 

As was shown in [8] , [5] polynomials of Rj type possesses not only orthogonality of the form 
(9) but also nice biorthogonality property with respect to some set of rational functions. 

In order to derive this biorthogonality propety in our case, we construct auxiliary polyno- 
mials 

T n (z) = P n+1 (z) -(z- a n+1 )P n (z). (11) 
Clearly, degree of these polynomials < n. More exactly, 

T n {z) = V n Z n + 0{z n - 1 ), 

where 

v n = a n+ i - a n H (12) 

«n+i a n 

In what folows we will assume that u n ^ 0. This means that degree of T n (z) is n and it is 
possible to introduce monic polynomials 

f n =T n (z)/Un = Z n + 0(z n - 1 ) (13) 

We have 

Theorem 1 Let P n (z) be Lagrange interpolation polynomials for the function F(z) and T n (z) 
defined by (11), (13). Define a set of rational functions 

v ^=a& <i4) 

Then Lagrange interpolation polynomials P n (z) and functions V n (z) form a biortogonal system 
in the folowing sense: 

Jr FjC) = n ( ' 

where the contour V should be chosen such that interpolation points ao, . . . lie inside the 
contour (N > max(n, m + 1)), and the function 1/F(z) is regular inside and on the contour. 

Proof. Assume first that m < n — 1. Then we have, obviously 

P n (()V m (0d( _ f Pn(()f m (()(z - a m+2 )(z - a m+3 ) ...(z- a n )d( _ 



r r n ((,)V m ((,)d(, = r 
Jr F(C) Jr 



F(C) Jr u; n+1 (z)F(() 

Pn(C)q n -i(z)dC 



/ r ^%P^o (1.) 



where q n -i(z) is a polynomial of degree < n — 1 and in the last equality in (16) we used 
orthogonality property (9). 
If m — n — 1 then 

p„(C)K-i(CR _ f P„(C)T»-i(C)dC _ 



ir F(C) ir 



*X0 ^ w^i(z)F(C) 



again by (9) 

If m > n then we can write down 



Pn{C)V m {QdC _ f P n (0(P m+1 (() - (( - a m+1 )P m (())d( 



Jr F(() ir 



/•'(0 ^ c m+2 (C)F(C) 

P n (C)A»+l(CK /" Pn(()Pm(()d( 



Ir c m+2 (C)F(C) ^ u m+1 (C)F(0 

because both terms in the last relation vanish due to (9) for n < m. 
Finaly, consider the case m = n. We have 

f Pn(QK(CR _ f Pn(0f n (()d( 

ir F(C) "ir u n+1 (0F(C) " 

Pn(()Pn + l(0d( f P*(()d( 







ir w n+2 (C)-F(C) ir Uni 



Ir u n+2 (C)F(0 Jru n+1 (()F(0 
Jr a n ^n(-2) 



The theorem is proven. 

The biorthogonality property can be rewritten in another form 



N 



V lim ((z - a s )P n (z)V m (z))/A s = 5 nm /a n , (17) 

_ Z Cta 



2 

s=0 



where N is any positive integer such that iV > max(n, m + 1). 

This result allows one to find coefficients in expansion of a given polynomial Q(z) of 
degree n in terms of interpolation polynomials P n (z) 

n 
fc=0 

Indeed, from (17) we have 

6 = {2m)- 1 a k [ Q(()V k (()d( (18) 



Consider an example. For the exponential function F(z) = exp(hz) (with an arbitrary 
nonzero real parameter h) choose uniform grid of the interpolation points a k — k, k = 0, 1, . . .. 
We then have (cf. [3]) 

Pn{z) = ± { -^{l-e h ) k (19) 

fc=0 K - 
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where (b) k — b(b + 1) . . . (b + k — 1) is the Pochhammer symbol. From (19) it is found 

a n = [e , iJ (20) 
n! 

Construct auxiliary polynomials T n {z) = P n+ \{z) — (z — a n+ i)P n (z). For the leading coefficient 
v n of the polynomials T n {z) = v n z n + 0(z n ~ l ) we have from (12) and (20) 



e h - 1 



7^0 



So polynomials T n (z) are indeed of degree n and for monic polynomials f n (z) = T n {z)/v n it is 
not difficult to obtain a rather attractive closed foormula 

T n {z)= ^ h + l) l 2Fl (- n >- Z ;l-e h ) (21) 
v ; (e h - l) n V-l -n ) K ' 

Thus for rational corresponding rational functions V n (z) we have from (14) 

Using standard transformation formulas for the Gauss hypergeometric function [1], we can 
present the functions V n (z) in a slightly different form 

Vn{z) = - — L - 2 F 1 (~ U, ~ Z ;e h ) (23) 

v ; (1 - e h ) n z(z - 1) V 2-z J K ' 

Thus rational functions V n (z) form a biorthogonal set with respect to the Lagrangean interpo- 
laiton polynomials (19): 

J r Pn(()V m (() eM-KW = S nm (24) 

where T is an arbitrary contour containing the points 0,1,..., max(n, m + 1) inside. 

Note finally, that recurrence relation (6) completely characterizes the Lagrange interpolation 
polynomials P n (z). More exactly, we have the 

Theorem 2 Assume that a set of monic n-th degree polynomials P n {z) satisfies recurrence 
relation (6) with initial conditions (7), where parameters a n , a n , n = 0, 1, . . . are arbitrary with 
the restrictions that all ai are distinct: a« ^ aj, for i ^ j and all a n are nonzero a n ^ 0, n = 
0, 1, . . .. Then polynomials P n (z) = a n P n (z) satisfy interpolation condition P n (a k ) = A k , k = 
0,1, ... ,n for all n = 0, 1, . . where 

n 

A n = P n (a n ) = J2 a s u s (a n ), n = 0, 1, 2, . . . (25) 

s=0 
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Proof. From recurrence relation (6) and initial conditions (7) it can be easily found 

P n+1 (z) —P n (z) = u n+1 (z), n = 0, 1, . . . (26) 

Hence for P n (z) = a n P n (z) we have the conditions 

P n+1 (a k ) = P n (a k ), k = 0, 1, . . . , n, n = 0, 1, 2, . . . 

From these relations, by induction, we obtain 

P n {a k ) = P k {a k ) = A k , k = 0, 1, . . . ,n 

Thus interpolation conditions are fulfilled. Expression (25) for A k follows then from Newton 
formula (3). 
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